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i!B ■ ABSTRACT 



Recent experimental data on r decays are used to reconstruct the difference in 
hadronic spectral densities with vector and axial-vector quantum numbers. The 
saturation of Das-Mathur-Okubo and Weinberg sum rules is studied. Two methods 
of improving convergence and decreasing errors are applied, and good agreement 
with the predictions of current algebra and chiral perturbation theory is observed. 
The resulting value of the pion polarisability is aE = (2.64 ± 0.36) 10~^ fm^. 
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1 Introduction 



The electric polarisability of the charged pion a e can be inferred from the amphtude 
for low energy Compton scattering 7 + 7r+ 7 + 7r+. This amplitude cannot be 
measured at low energies directly, but can be determined from measurements on 
related processes like ttN —>■ ttN'^, jN — * '-/Nn and 77 —>■ txtx. The measured values 
for aE, in units of 10"^ fm^ are 6.8 ± 1.4 ± 1.2 H, 20 ± 12 and 2.2 ± 1.6 0, 
respectively. 

Alternatively, the polarisability can be predicted theoretically by relating it to 
other quantities which are better known experimentally. In chiral perturbation 
theory, it can be shown that the pion polarisability is given by 

a Fa , . 

where F-^ is the pion decay constant and Fa is the axial structure-dependent form 
factor for radiative charged pion decay n — > ez/7 |^ . The latter is often re-expressed 
in the form 

Fa = jFv, (2) 

because the ratio 7 can be measured in radiative pion decay experiments more accu- 
rately than Fa itself, while the corresponding vector form factor Fy is determined to 
be 10 Fv = 0.0131 ± 0.0003 by using the conserved vector current (CVC) hypoth- 
esis to relate vr ewy and 7r° 77 decays. 7 has been measured in three vr ei/j 
experiments, giving the values: 0.25 ± 0.12 0; 0.52 ± 0.06 and 0.41 ± 0.23 [|. 
The weighted average 7 = 0.46 ± 0.06 can be combined with the above equations to 
give 

aE = (2.80 ± 0.36) 10"^ fm^ . (3) 

This result is often referred to as the chiral theory prediction for the pion polar- 
isability [Q. However a^;, or equivalently Fa, can be determined in other ways. In 



particular, the latter occurs in the Das-Mathur-Okubo (DM0) sum rule 

I = Fr-^-FA, (4) 



^Oiir definitions of Fy and Fa diS^er from those used by the Particle Data Group g by a factor 
of two. 
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where 

/ = / ^Pv-a{s) , (5) 

with pv^a{s) = pv{s) — pa{s) being the difference in the spectral functions of the 
vector and axial-vector isovector current correlators, while (r^) is the pion mean- 
square charge radius. Using its standard value (r^) = 0.439 ± 0.008 fm^ [0] and 
eqs. i^, (0) one gets: 

/mfo = (26.6±1.0)-10^=^ (6) 

Alternatively, if the integral / is known, eq. (^) can be rewritten in the form of 
a prediction for the polarisability: 

a ({vD I 



Recent attempts to analyse this relation have resulted in some contradiction with 
the chiral prediction. Lavelle et al. |l2| use related QCD sum rules to estimate the 
integral / and obtain ue = (5.60 ± 0.50) 10~^ fm^. Benmerrouche et. al. [jl3| apply 
certain sum rule inequalities to obtain a lower bound on the polarisability (|^) as a 
function of (r^). Their analysis also tends to prefer larger oe and/or smaller (r^) 
values. 

In the following we use available experimental data to reconstruct the hadronic 
spectral function pv-a{,s), in order to calculate the integral 

/■So ds 

/o(so) = / —Pv-a{s) . (8) 

for So — M^, test the saturation of the DM0 sum rule (^ and its compatibility with 
the chiral prediction d). We also test the saturation of the first Weinberg sum rule 



0: 

W,{so) ^ / dspv-A{s) ; W,{so) L„^^= (9) 

and use the latter to improve convergence and obtain a more accurate estimate for 
the integral /: 

r^o ds 



f^o (is 

h{so) = / —Pv-a{s) 

JAml S 



So 



F^- dspv-Ais) 

J Ami 



(10) 
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Here the parameter (3 is arbitrary and can be chosen to minimize the estimated error 
in/i m- 

Yet another way of reducing the uncertainty in our estimate of / is to use the 



Laplace-transformed version of the DM0 sum rule |17 

I^[M^) = Fl^^ - Fa 



'ID 



with being the Borel parameter in the integral 



f ds f —s 

C,{0^) Cs{0') 



exp I — ) Pv-a[s) + 



F^ 
M2 



(12) 



6M6 24M8 

Here Cq{0^) and Cs{0^) are the four-quark vacuum condensates of dimension 6 and 
8, whose values could be estimated theoretically or taken from previous analyses 



15, 18 



All the three integrals (|]), ([T0| ) and (|T2D obviously reduce to (g) as Sq, M — >• oo. 



2 Evaluation of the spectral densities 

Recently ALEPH published a comprehensive and consistent set of r branching frac- 
tions [|l^], where in many cases the errors are smaller than previous world averages. 
We have used these values to normalize the contributions of specific hadronic final 
states, while various available experimental data have been used to determine the 
shapes of these contributions. Unless stated otherwise, each shape was fitted with 
a single relativistic Breit-Wigner distribution with appropriately chosen threshold 
behaviour. 

2.1 Vector current contributions. 



A recent comparative study ||2^ of corresponding final states in r decays and e"*"e~ 
annihilation has found no significant violation of CVC or isospin symmetry. In order 
to determine the shapes of the hadronic spectra, we have used mostly r decay data, 
complemented by e~^e~ data in some cases. 

7r~7r° : BR = 25.30 ± 0.20% []Ty[, and the s-dependence was described by the 
three interfering resonances p(770), p(1450) and /9(1700), with the parameters taken 
from IBi and \Ul . 



Stt^tt^ : BR = 4.50 ± 0.12%, including tt uj final state |T^. The shape was 



determined by fitting the spectrum measured by ARGUS p2 |. 

7r~37r° : BR = 1.17 ± 0.14% [19|. The s-dependence is related to that of the 
reaction e"*"e~ 27r+27r^. We have fitted the latter measured by OLYA and DM2 



2|. 



Gtt : various charge contributions give the overall BR = 0.13 ± 0.06% [1^, fairly 
close to CVC expectations . 



TT ri] : BR = 0.17 ± 0.03% The s-dependence was determined by fitting 



the distribution measured by CLEO |23 



K'K^ : BR = 0.26 ± 0.09% Again, the fit of the CLEO measurement |25 
was performed. 

2.2 Axial current contributions. 

The final states with odd number of pions contribute to the axial- vector current. 
Here, r decay is the only source of precise information. 

TT^ : BR = 11.06 ± 0.18% ||T9|. The single pion contribution has a trivial s- 
dependence and hence is explicitly taken into account in theoretical formulae. The 
quoted branching ratio corresponds to F^^ = 93.2 MeV. 

37r± and 7r-27r° : BR = 8.90 ± 0.20% and BR = 9.21 ± 0.17%, respectively JIU . 
Theoretical models assume that these two modes are identical in both shape and 
normalization. The s-dependence has been analyzed in where the parameters 
of two theoretical models describing this decay have been determined. We have used 
the average of these two distributions, with their difference taken as an estimate of 
the shape uncertainty. 

37r±27r° : BR = 0.50 ± 0.09%, including h'tt^uj final state 111. The shape was 



fitted using the CLEO measurement |28 

57r± and tt^Atv'^ : BR = 0.08 ± 0.02% and BR = 0.11 ± 0.10%, respectively p 



We have assumed that these two terms have the same s-dependence measured in 
2|. 
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2.3 KKtt modes. 

KK-K modes can contribute to both vector and axial-vector currents, and various 
theoretical models cover the widest possible range of predictions . 



According to |1T9|, all three KK-k modes {K K^-r- , R-R^tt^ and R-R+n-) add 
up to BR.{RRtt) = 0.56 ±0.18%, in agreement with other measurements (see p5[|). 
The measured s-dependence suggests that these final states are dominated by R*R 
decays We have fitted the latter, taking into account the fact that due to parity 
constraints, vector and axial- vector R*R terms have different threshold behaviour. 
A parameter ^ was defined as the portion of RRn final state with axial-vector 
quantum numbers, so that 

BR(RRtt)v = (1-0 BR(RRn) 

BR(RRtt)a = ^BR(RR7t). (13) 

3 Results and conclusions 

The resulting spectral function is shown in Fig.l. The results of its integration 
according to (^ are presented in fig. 2 as a function of the upper bound sq. One can 
see that as sq increases, Jq converges towards an asymptotic value which we estimate 
to bei 

/o = /o(oo) = (27.5 ± 1.4 ± 1.2) ■ 10"^ (14) 

in good agreement with the chiral value (|^). 

The saturation of the Weinberg sum rule (|^) is shown in fig. 3. One sees that 
the expected value is well within the errors, and ^ ~ 0.25 ^ 0.3 seems to be 
preferred. No significant deviation from this sum rule is expected theoretically |]l6l , 



so we use (0) to calculate our second estimate of the integral /. The results of this 
integration are presented in fig. 4, with the asymptotic value 

h = h{oo) = (27.0 ± 0.5 ± 0.1) ■ 10-^ (15) 

corresponding to (3 ~ 1.18. One sees that the convergence has improved, the errors 
are indeed much smaller, and the (^-dependence is very weak. 



^In the following, the first error corresponds to the quadratic sum of the errors in the branching 
ratios and the assumed shapes, while the second one arises from to the variation of ^ in the interval 
0.5 ±0.5. 
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Figure 1: Difference of vector and axial-vector hadronic spectral densities. In figs. 1-5: the three 
curves correspond to ^ = 0, 0.5 and 1 from top to bottom; the errors originating from the shape 
variation and those coming from the errors in the branching fractions are roughly equal and have 
been added in quadrature to form the error bars, shown only for ^ = 0.5. 



Now we use to obtain our third estimate of the spectral integral. The integra- 
tion results are plotted against the Borel parameter in fig.5, assuming standard 
values for dimension 6 and 8 condensates. The results are independent of for 
> 1 GeV^, indicating that higher order terms are negligible in this region, and 
giving 

h = hioo) = (27.2 ± 0.4 ± 0.2 ± 0.3) 10"^ , (16) 

where the last error reflects the sensitivity of (^ to the variation of the condensate 
values. 

One sees that these three numbers (|1^ - (0) are in good agreement with each 
other and with the chiral prediction (^). Substitution of our most precise result (|T5|) 
into (0) yields for the standard value of the pion charge radius quoted above: 



oe = (2.64 ±0.36) 10"^ fm^ 



(17) 



in good agreement with (|^) and the smallest of the measured values, 0. Note that 
by substituting a larger value (r^) = 0.463 ± 0.006 fm^ one obtains aE = 
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Figure 2: Saturation of the DM0 sum rule integral (8). The thick dashed line is the chiral 
prediction for the asymptotic value (6) and the dotted lines show its errors. 
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Figure 3: Saturation of the first Weinberg sum rule (9). The dashed line shows the expected 
asymptotic value F?^. 
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Figure 4: Saturation of the modified DM0 sum rule integral (10). The chiral prediction is also 
shown as in fig. 2. 
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Figure 5: The Laplace-transformed sum rule (12) as a function of the Borel parameter iV/^, 
compared to the chiral prediction. 
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(3.44 ± 0.30) 10~^ fm^, about two standard deviations higher than (^. 

In conclusion, we have used recent precise data to reconstruct the difference in 
vector and axial-vector hadronic spectral densities and to study the saturation of 
Das-Mathur-Okubo and the first Weinberg sum rules. Two methods of improving 
convergence and decreasing the errors have been used. Within the present level of 
accuracy, we have found perfect consistence between r decay data, chiral and QCD 
sum rules, the standard value of (r^), the average value of 7 and the chiral prediction 
for aE- 

Helpful discussions and correspondence with R. Alemany, R. Barlow, M.Lavelle 
and P. Poffenberger are gratefully acknowledged. 
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